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SYMBOLS  AND  ABBREVIATIONS 


Notation 


Definition 


E 


A  Pc 
r 

u 

V 

Vr  and  Vg 

Y 
z 


£ip,  Szp,  and 


Ap 

APo,  and 


p>p  and  pt 


Oh 


internal  energy 

pressure,  compaction  pressure 
radial  coordinate 
penetration  rate 

velocity  of  the  rigid  portion  of  the  penetrator 
radial  and  axial  compnent  of  the  velocity 
von  Mises  yield  stress 
axial  coordinate 

principle  plastic  strain  components 

size  of  the  plastic  zone  in  front  of  the  penetrator 

density,  initial  density,  and  reference  density, 
respectively 

density  of  the  penetrator  and  the  target, 
respectively 

normal  stress  at  a  boundary 
porosity,  initial  porosity 


INTRODUCTION 


When  a  high-velocity  projectile  impacts  a  target,  the  penetration  is  accomplished  by 
displacing  the  target  material  radially,  and  as  a  result,  a  tunnel  shaped  crater  is  formed  in 
the  target.  The  depth  and  the  diameter  of  the  produced  hole  depend  upon  the  projectile 
velocity  and  the  material  properties  of  the  projectile  and  the  target.  In  particular,  the  "flow" 
of  the  target  material  under  the  pressure  exerted  by  the  projectile  is  very  sensitive  to  the 
nature  and  mode  of  failure  of  the  target  material.  Concrete  is  a  very  complex  material, 
consisting  of  mineral  aggregate  bound  by  cement  paste  containing  a  large  amount  of  water 
and  voids.  As  a  result,  it  is  only  recently  that  the  behavior  of  this  complex  material  in 
response  to  high-velocity  penetration  has  begun  to  be  understood. 

Only  a  limited  number  of  three-dimensional  concrete  penetration  calculations  have 
been  reported  in  the  literature.  These  calculations  (refs  1  to  3)  pertain  to  simulations  of  low 
velocity  impacts  and  employ  codes  based  on  a  Lagrangian  formulation  of  the  governing 
equations.  Since  in  case  of  low  velocity  impacts,  the  projectile  is  subjected  to  relatively 
small  deformations  and  remains  relatively  rigid  and  nondeformed,  the  Lagrangian  approach 
can  offer  some  computational  advantages.  However,  even  when  the  deformations  of  the 
penetrator  are  relatively  small,  the  target  elements  adjacent  to  the  penetrator  front  are 
deformed  severely  and  cannot  "flow"  around  the  penetrator.  To  overcome  this  difficulty 
requires  frequent  re-zoning  of  the  target  material  in  the  neighborhood  of  the  penetrator 
front.  The  re-zoning  adds  to  the  computation  time  and  allows  for  variability  in  the  results, 
depending  on  the  type  of  the  re-zoning  algorithm  and  the  frequency  of  re-zoning  (ref  1 ).  A 
number  of  issues  and  techniques  applicable  to  the  penetration  calculations  with  Lagrangian 
codes  are  discussed  in  reference  1 . 

The  concrete  constitutive  models  employed  in  all  the  above  mentioned  studies  exhibit 
limited  features  such  as  pressure-dependent  yield  surfaces  (refs  1  and  2),  pressure  and 
stain-rate-dependent  yield  surfaces,  and  pressure-dependent  bulk  and  shear  moduli  (ref  3). 
Hydrocode  implementation  of  the  existing  advanced  concrete  models  encounter  a  number 
of  computer  code  development  challenges.  This  is  mainly  due  to  the  extensive 
computational  resources  required  for  the  calculations.  Recently,  a  number  of  advanced 
models  that  are  capable  of  adequately  modeling  the  complex  dynamic  behavior  of  concrete 
has  been  reported  in  the  literature.  Among  these  models  are  the  Bicanic  and  Zienkiewicz 
model  (refs  4  and  5)  and  the  porous  material  models  (refs  6  and  7). 

In  this  work,  state-of-the-art  advanced  models  for  concrete  are  applied  to  a  problem  of 
high-velocity  penetration,  enabling  the  prediction  of  the  target  cratering  and  hole  shape.  A 
schematic  of  a  crater  produced  in  a  concrete  target  (ref  8)  is  shown  in  figure  1 .  The 
entrance  portion  of  the  crater  is  produced  by  the  spallation  of  concrete  in  the  vicinity  of  the 
impact  site.  The  crater  entrance  is  typically  wide  and  shallow,  rapidly  evolving  into  a  well 
rounded  and  slightly  tapered  tunnel.  The  penetration  prediction  capability  of  various 
constitutive  models  for  concrete,  which  were  studied  in  this  work,  was  assessed  using  the 
shape  of  the  latter  portion  of  the  crater  only.  The  accuracy  of  the  penetration  prediction 


1 


capability  of  these  models  was  established  by  comparing  the  numerical  results  to 
experimental  data  (ref  8). 

When  the  rate  of  penetration  in  concrete  is  below  the  sonic  velocity,  a  zone  of 
comminuted  concrete  is  formed  in  front  of  the  penetrator,  and  the  extent  of  this  zone  is 
influenced  by  the  dynamic  properties  of  the  target  material.  A  number  of  experimental 
techniques  for  acquisition  of  the  subsonic  target  material  flow  in  front  of  the  penetrator  are 
reported  in  the  literature  (refs  9  and  10).  Unique  data  of  the  dynamic  behavior  of  concrete 
have  been  obtained  and  reported  earlier  (ref  8)  using  a  sequence  of  break  gages.  These 
data,  as  well  as,  numerous  experimental  and  analytical  studies  of  penetration  with  high 
length  over  diameter  (LA))  ratio  projectiles  indicate  that  the  trajectory  of  the  front  end  of  the 
penetrator  consists  of  two  distinct  regions:  a  long  region  of  approximately  steady  state 
penetration  and  a  short  region  of  deceleration.  A  photograph  of  a  sectioned  target  from 
these  experiments  in  shown  in  figure  2.  Break  gages  that  were  originally  plain  are 
uniformly  bulged  out  in  the  direction  of  the  projectile  motion.  Note  that  all  gages  are 
approximately  parallel  to  each  other,  which  may  lead  to  a  conjecture  regarding  the  steady- 
state  motion  of  the  penetrator  in  the  target. 

The  technique  for  determining  the  extent  of  the  plastic  zone  in  front  of  the  projectile  is 
based  on  the  fact  that  the  projectile  penetration  rate  in  the  steady-state  region  depends 
strongly  on  the  densities  of  the  projectile  and  the  target,  and  only  weakly  on  their  strengths. 
This  allows  an  accurate  prediction  of  this  region  of  the  trajectory  without  detailed 
knowledge  of  the  constitutive  behavior  of  the  target  material.  In  the  experiments  by  Gold 
et.  al.  (ref  8),  break  gages  were  placed  along  the  projectile  trajectory  in  a  region  of  steady- 
state  penetration.  Voltages  applied  to  the  gages  were  monitored  for  sudden  interruptions 
(i.e.,  for  'breaks"  of  the  gages),  which  were  associated  with  the  time  of  arrival  of  the  front  of 
the  zone  of  fractured  concrete.  Once  the  trajectory  of  the  projectile  front  end  is  calculated, 
the  extent  to  the  zone  of  comminuted  concrete  can  be  measured  as  the  difference  between 
the  trajectory  and  the  regression  of  trigger  times  of  the  break  gages. 

In  the  present  work,  the  numerical  simulations  have  concentrated  on  modeling  an 
earlier  experiment  involving  a  19-cm  long  and  1.3-cm  dia  copper  projectile  impacting  a 
plain  concrete  target  with  velocity  Vo  =  0.1 836  cm/jus  (ref  8).  The  projectile  was  carefully 
aimed  along  the  target  axis,  and  x-ray  records  confirmed  a  90-deg  incident  angle  at  the 
impact.  These  experimental  conditions  simulated  an  axially  symmetrical  configuration,  and 
the  analytical  study  was  reduced  to  a  three-dimensional  axially  symmetrical  problem.  This 
work  presents  results  of  numerical  calculations  using  various  constitutive  models  for 
concrete. 


NUMERICAL  SCHEME  AND  MATERIAL  MODELS 

The  numerical  simulations  presented  in  this  work  were  performed  using  the  CALE  (ref 
1 1)  computer  program.  The  CALE  is  the  plane  two-dimensional  and  three-dimensional 
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axisymmetric  hydrodynamics  code  based  on  an  Arbitrary  Lagrangian  Eulerian  (ALE) 
formulation  of  the  governing  equations.  It  is  capable  of  modeling  multi-material  flows  and 
allows  for  discontinuous  velocities  at  material  interfaces.  In  addition  to  the  specification  of 
the  problem  geometry  and  initial  and  boundary  conditions,  the  equations  of  state  (EOS) 
and  the  constitutive  equations  for  the  target  and  the  projectile  have  to  be  specified  before 
the  solution  procedure  can  be  initiated. 

The  constitutive  equation  for  the  projectile  material  adopted  here  is  that  of  Steinberg- 
Guinan  (ref  12),  which  specifies  the  shear  modulus  and  the  yield-strength  as  a  function  of 
pressure,  temperature,  and  equivalent  plastic  strain.  The  EOS  that  was  used  to  model  the 
hydrodynamic  response  of  copper  was  the  standard  linear  polynomial  approximation 
usually  employed  for  metals  (ref  11).  Since  the  copper  projectiles  were  manufactured  from 
a  standard  stock  of  oxygen-free  high  conductivity  copper  (OFHC),  the  analysis  relied  on  a 
standard  set  of  parameters  for  this  material  from  reference  13.  Details  of  these  models  and 
their  implementation  can  be  found  in  reference  11. 

Two  constitutive  models  were  used  to  represent  the  behavior  of  concrete: 

•  An  elastic-plastic  model  with  von  Mises  yield  surface 

•  An  elastic-plastic  model  with  pressure-dependent  yield  surface 

The  sensitivity  of  the  impact  response  of  concrete  to  the  selection  of  an  EOS,  a  constitutive 
model,  and  relevant  parameters  was  studied  in  this  work.  Details  regarding  these  models 
and  results  of  the  calculations  are  discussed  later. 


NUMERICAL  MODEL 

The  geometry  and  boundary  conditions  are  shown  in  figure  3.  The  domain 
considered  is  axisymmetric  in  r-z  space.  Domain  T^T2T3T4  represents  the  target,  P1P2P3T1 
is  the  spherical-nose  projectile,  and  the  rest  of  the  domain  E^E2T2T\  is  filled  with  air.  Eight 
0.03-cm  thick  break  gages  were  initially  modeled  using  three  computations  cells,  each 
0.01 -cm  thick,  filled  with  copper.  After  some  initial  numerical  experimentation  it  became 
clear  that  these  0.01 -cm  cells  were  resulting  in  very  small  time  steps,  which  controlled  the 
overall  speed  of  calculations.  Therefore,  in  order  to  decrease  the  otherwise  prohibitively 
excessive  computational  time,  the  number  of  cells  across  the  thickness  of  the  break  gages 
was  reduced  to  one. 

The  following  boundary  conditions  were  applied  along  the  boundary  of  domain 

E1E2T3T4: 


Vr  =  V2=0  along  E1E2 
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Oh  =  0  along  E2T3  and  T3T4 

Vr  =  0  along  E^  T4  (line  of  axial  symmetry) 

The  portion  of  the  mesh  A-[A2TzT4  was  run  in  ALE  mode,  but  the  nodes  within  the 
domain  E1E2A2A1  were  reset  to  their  original  position  after  each  cycle  of  calculations. 
Attempts  to  run  the  domain  E1E2A2A1  in  ALE  mode  resulted  in  severe  deformations  and 
tangling  up  of  the  ALE  mesh  due  to  the  debris  of  concrete  ejected  from  the  crater  and 
moving  away  from  the  target  surface.  This  significantly  reduced  the  time  step,  slowed  down 
calculations,  and  often  caused  termination  of  the  run.  The  present  work  primarily 
concentrated  on  simulating  the  formation  of  the  crater,  rather  than  on  a  study  of  a  spallation 
off  the  front  surface  of  the  target,  which  would  require  a  detailed  analysis  of  the  motion  of 
this  debris.  Therefore,  spacing  the  nodes  in  domain  E1E2A2A1  coarsely  and  "running"  them 
in  the  Eulerian  mode  significantly  sped  up  the  calculations.  Thus,  a  significant  increase  in 
the  computational  speed  was  achieved,  although,  at  an  acceptable  loss  of  accuracy  of  the 
analysis  in  that  region.  In  a  few  trial  runs,  the  portion  of  the  grid  contained  in  E1E2A2E2  was 
omitted,  the  free  boundary  was  set  along  R2A2  and  the  boundary  conditions  w  =  Vz  =  0  were 
imposed  along  R^R2.  This  allowed  a  slight  reduction  in  the  number  of  computational  zones 
and  sped  up  the  calculations.  However,  it  was  determined  that  the  grid  motion  algorithms 
were  very  sensitive  to  variations  in  the  values  of  the  parameters  appearing  in  the  material 
models,  which  would  in  turn  complicate  the  course  of  the  parametric  studies.  Thus,  this 
approach  was  abandoned. 

A  non-uniform  distribution  of  nodes  in  both  the  radial  and  axial  directions  was 
employed  in  the  discretization  model.  The  computational  mesh  for  all  the  zones  was 
designed  based  on  the  rationale  that  the  best  accuracy  of  analysis  would  be  achieved  with 
fine  degrees  of  discretization  in  the  areas  subjected  to  high  gradients,  and  a  coarser  mesh 
in  regions  with  moderate  gradients.  The  finest  radial  grid  size  Armm  was  set  along  the  path 
of  motion  of  the  projectile  where  the  material  deformations  are  the  most  severe.  The  radial 
grid  size  gradually  increased  toward  the  maximum  Aw  at  the  lateral  boundary.  In  the 
axial  direction,  zones  1  through  4  (fig.  3)  denote  regions  with  different  axial  grid  sizes  Az. 
Nodes  in  zone  2  were  spaced  uniformly.  In  zones  1 ,  3,  and  4  axial  grid  sizes  increased 
gradually  from  the  smallest  AZmin  at  one  boundary  to  the  largest  Azmax  at  the  opposite 
boundary.  The  axial  size  of  zone  3  was  selected  according  to  the  expected  depth  of 
penetration.  In  zone  4,  which  is  beyond  the  expected  depth  of  penetration,  nodes  were 
spaced  coarsely  (in  the  axial  direction).  In  order  to  resolve  the  formation  of  the  crater 
entrance,  the  grid  in  zone  1  was  finely  meshed  close  to  the  target  surface,  while  the  size  of 
the  computational  mesh  gradually  increased  toward  the  boundary  E1E2. 

The  effect  of  the  grid  size  on  the  overall  accuracy  of  the  calculations  is  summarized  in 
table  1 .  In  all  cases  the  computational  mesh  had  a  similar  configuration,  but  the  total 
number  of  the  computational  elements  varied  for  29,895  to  9,437.  Reduction  of  the  total 
number  of  elements  from  29,895  to  10,971  was  achieved  by  increasing  the  zone  of  the 
elements  in  zones  2  and  3  by  a  factor  of  two,  which  cut  the  required  CPU  time  of  SUN 
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4/490  for  a  sample  200//S  long  event  from  1 02.2  to  1 8.4  hrs.  Reduction  of  the  number  of 
elements  from  10,971  to  9,437  was  achieved  by  changing  the  number  of  cells  across  the 
thickness  of  the  gages  from  three  to  one,  which  increased  the  computational  speed  by  a 
factor  of  about  two.  Since  this  study  concentrated  on  the  concrete,  this  approach  had  no 
appreciable  effect  on  the  accuracy  of  the  calculations.  Changes  in  the  modes  for 
controlling  grid  motion  from  mixed  Euler/ALE  to  purely  Eulerian  analysis  significantly 
affected  the  speed  of  calculation  reducing  if  (for  the  Eulerian  analysis)  by  a  factor  of  about 
three.  Calculated  crater  profiles  compared  at  200^s  were  nearly  the  same  for  all  cases 
with  only  slight  differences  in  the  resolution  of  the  crater  lip  formed  by  the  ejected  debris  of 
concrete.  Since  the  accuracy  of  analysis  was  not  affected  either  by  varying  the  number  of 
computational  cells  or  by  changing  modes  of  grid  motions,  the  mesh  configuration  resulting 
in  the  lowest  CPU  time  (pci0_02)  was  adopted  for  all  parametric  studies  presented  in  this 
work. 


RESULTS 


Concrete  Targets:  The  EOS  Model 

An  estimate  for  the  low  bound  of  pressures  exerted  on  the  concrete  at  the  front  end  of 
the  projectile  impacting  with  velocity  Vo  =  0.1836  cm/^iS  can  be  obtained  from  the  ideal  one¬ 
dimensional  hydrodynamic  theory  of  penetration  (ref  13).  Indeed,  the  following  equation; 


P  =  Tpp{v  -  uY  =  Yptu^  (1) 

with  pp  =  8.9glcm^  for  copper  projectile  and  =2.24g/cm^  for  concrete,  yields 
u  =  0.\22cm  t  ps ,  which  gives  a  value  for  the  pressure  of  p  =  1 6.7  Kbar.  More  realistic 
estimates  of  the  pressure  are  obtained  from  three-dimensional  axisymmetrical  elastic- 
plastic  analysis  which  yields  values  of  the  order  of  20  to  27  Kbar  (fig.  4).  At  these  high 
pressures,  which  exceed  the  static  unconfined  strength  of  concrete  by  almost  two  orders  of 
magnitude,  the  effects  of  compressibility  of  the  target  material  are  of  significant  importance. 

Experimental  shock-wave  studies  for  concrete  conducted  by  Gregson  (ref  14)  reveal 
significant  non-linearity  in  Hugoniot  data  for  particle  velocities  below  0.1  cm/ps.  Non¬ 
linearity  in  the  Hugoniot  can  be  attributed  to  physical  mechanisms  of  micro-cracking  and 
pore  closing.  Read  and  Maiden  (ref  15),  using  a  porous  EOS  model,  were  able  to 
reproduce  the  essential  features  of  the  experimental  shock  velocity  vs  particle  (material) 
velocity  data.  The  particle  velocities  experienced  by  the  concrete  in  the  earlier  penetration 
experiments  (ref  8)  are  in  the  very  same  range.  Therefore,  applying  the  porous  EOS  model 
is  appropriate  for  describing  the  constitutive  response  of  concrete. 
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Following  Tipton  (ref  6),  the  porous  EOS  model  assumes  that,  while  in  compression, 
the  porosity  is  not  changed  until  the  pressure  exceeds  the  compaction  pressure,  Pc.  In 
tension,  the  porosity  is  not  increased  until  the  pressure  drops  below  -  Pc.  The  compaction 
pressure  is  assumed  to  be  a  function  of  the  porosity  ^and  has  the  following  form: 


pc  =  po\og{(p)  (2) 

The  total  pressure  of  the  porous  material  is  given  by 

p  =  {\-(p)ps(,ps)  (3) 

where  <p  is  the  porosity,  and  Ps  and  ps  are  the  pressure  and  the  density  of  the  solid  phase, 
respectively.  The  EOS  for  the  solid  phase  is  taken  in  a  polynomial  form: 

p- Ao-\- A\p-\- Aip  +  Azp^ +{Bo  +  B\p  +  Bip^)E  (4) 

where  p  =  pt pr-\,pr  =  pol{\-(pb)  is  reference  density,  and  is  the  initial  porosity. 

Further  details  on  this  model  and  its  implementation  in  GALE  can  be  found  in  reference  11. 

A  comparison  of  various  flow  field  variables  calculated  with  solid  (i.e.,  directly  via 
equation  4  and  porous  (i.e.,  using  equations  3  and  4)  EOS  models  are  presented  in  figures 
4  to  1 1 .  The  values  of  parameters  for  these  constitutive  models  are  given  in  table  2.  The 
coefficients  appearing  in  the  polynomial  EOS  form  were  adopted  from  reference  1 1  and  are 
based  upon  the  Hugoniot  data  from  reference  14.  In  all  cases,  the  calculations  were 
carried  out  up  to  SOOps.  All  calculations  show  that  the  projectile  almost  completely  erodes 
at  300  to  400ps,  and  the  penetration  process  essentially  ceases  at  approximately  SOOps. 
Figure  5  compares  the  experimental  data  with  the  calculated  final  hole  profiles  at  500ps. 
Regardless  of  the  EOS  model  (as  well  as  other  models  to  be  considered  later),  the 
entrance  portion  of  the  crater  is  predicted  poorly.  All  numerical  calculations  presented  in 
this  work  failed  to  reproduce  not  only  the  size  of  the  entrance,  but  mostly  importantly,  even 
the  general  shape  of  this  portion  of  the  crater.  Since  the  entrance  portion  of  the  crater  is 
produced  by  spallation  of  concrete  in  the  vicinity  of  the  impact  site,  its  size  and  shape  are 
strongly  influenced  by  brittle  fracture,  a  phenomenon  that  was  not  modeled  in  the  present 
work.  Experiments  show  that  the  crater  entrance  is  typically  wide  and  shallow,  rapidly 
evolving  into  a  well  rounded  and  slightly  tapered  tunnel.  The  penetration  prediction 
capability  of  various  constitutive  modes  for  concrete  studied,  was  assessed  using  the 
shape  of  the  latter  portion  of  the  crater  only.  The  accuracy  of  the  penetration  prediction 
capability  of  these  models  was  established  by  comparing  the  numerical  results  to 
experimental  data  (ref  8).  Both  calculations  (fig.  5)  significantly  disagree  with  the 
experimental  data.  Formation  of  the  tunnel  portion  of  the  crater  is  accomplished  through 
plastic  flow  of  the  target  material,  driven  by  the  pressure  exerted  at  the  interface  between 
the  projectile  and  the  target.  Under  these  conditions,  this  material  is  always  under 
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compression  and  is  far  removed  from  the  region  of  the  crater  entrance  where  the  material  is 
in  tension  and  is  subjected  to  brittle  fracture.  Therefore  the  influence  of  the  brittle 
properties  of  concrete  on  the  formation  of  the  tunnel  portion  of  the  crater  is  negligible. 
Experimental  data  indicate  that  the  shape  of  this  portion  of  the  crater  is  similar  to  the  shape 
of  boreholes  in  ductile  materials,  such  as  metals,  etc.  (ref  8). 

Comparisons  of  results  obtained  with  the  solid  and  porous  EOS  models  reveal  that 
the  hole  profile  computed  with  the  porous  model  is  slightly  deeper  and  almost  twice  as  wide 
as  the  profile  predicted  with  the  solid  model.  From  the  point  of  view  of  a  one-dimensional 
modified  hydrodynamic  theory  of  penetration,  the  penetration  depth  is  a  function  of  the 
pressure  at  the  front  end  of  the  projectile.  Thus,  the  average  pressures  computed  with  the 
porous  EOD  model  were  slightly  lower  than  with  the  solid  EOS  model,  and  this  is  in 
agreement  with  the  Hugoniot  data  from  shock-wave  studies. 

Figure  4  presents  plots  of  isolines  of  pressures  computed  at  approximately  35jus. 
Although  most  of  the  penetration  process  is  quasi-steady  state,  the  calculated  pressures 
reflect  specific  wave  propagation  patterns  which  are  very  sensitive  to  material  models. 

Both  plots  shown  in  figure  4  are  taken  at  approximately  the  same  time;  however,  the 
computed  wave  propagation  patterns  reflect  slightly  different  wave  propagation  states. 

This  explains  the  fact  that  the  maximum  pressure  recorded  in  the  plot  for  the  porous  EOS 
model  is  higher  than  that  for  the  solid  EOS  model.  Nevertheless,  for  the  solid  EOS  model, 
the  overall  pattern  of  widely  spread  isolines  of  pressure  is  distinctly  different  from  pressure 
distributions  computed  with  the  porous  EOS  model.  This  essentially  confirms  the 
qualitative  features  that  are  expected  from  "snow  plow"  type  of  EOS  model. 

The  calculated  velocity  fields  are  less  prone  to  fluctuations  (due  to  stress  wave 
propagation  )  and  are  quasi-steady  state.  Figure  6  shows  profiles  of  the  axial  component 
of  velocity  along  the  centerline  at  time  35jus.  The  interface  between  projectile  and  concrete 
is  located  approximately  at  z  =  4.4  cm,  and  although  in  the  projectile  (z<z*)  the  velocity 
profile  is  the  same,  the  distribution  of  the  velocities  in  the  concrete  (z>z*)  is  quite  different. 
The  velocity  profile  calculated  with  the  solid  EOS  model  is  smooth  and  the  flow  spreads  out 
for  significant  distances  in  front  of  the  projectile.  The  situation  is  quite  opposite  for  the 
porous  EOS  model.  The  flow  is  contained  in  the  immediate  vicinity  of  the  penetrator  front 
and  is  rapidly  diminished  away  from  it.  The  same  fact  is  established  from  comparison  of 
the  radial  flow  fields  shown  in  figures  6  and  7.  As  with  the  axial  velocity  profiles,  for  the 
solid  EOS  model,  the  profile  of  the  radial  component  of  velocity  is  smooth  and  spread  out, 
while  for  the  porous  EOS,  the  flow  field  is  more  intense  and  confined.  Since  final  hole 
diameters  are  proportional  to  the  magnitude  of  the  radial  flow,  the  porous  EOS  model  is 
expected  to  result  in  significantly  wider  holes,  which  in  turn,  correlates  with  the  results 
presented  in  figure  5. 

Results  presented  in  figures  8  to  10  summarize  the  effects  of  the  EOS  model  on  the 
flow  field  structure  in  front  of  the  projectile.  Calculated  patterns  of  isolines  of  the  equivalent 
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plastic  strain  are  shown  in  figure  8.  The  equivalent  plastic  strain  is  defined  as: 


y  ~  3  ^2p)  i^2p  ^3p)  ^\p)  (5) 

where  8ip,  E2p,  and  esp  are  the  principal  plastic  strain  components,  and  reflects  the  extent  of 
elastic-plastic  state.  Isolines  of  a  small  y  mark  the  position  of  the  elastic-plastic  boundary  in 
concrete.  A  Comparison  between  these  two  plots  yields  that  the  sizes  of  plastic  zones  in 
front  of  the  projectile  are  remarkably  different  for  the  two  EOS  models.  A  plot  of  equivalent 
plastic  strain  along  the  streamline  r  =  0  is  shown  in  figure  9.  Consistent  with  the  previous 
results  for  the  porous  EOS  model,  the  plastic  flow  is  more  intense  and  confided  to  a  short 
zone  in  front  of  the  projectile.  The  position  of  the  elastic-plastic  boundary  is  distinct  and  the 
approximate  size  of  the  plastic  zone  in  concrete  in  front  of  the  projectile  can  be  measured 
directly  from  the  plot,  yielding  Ap » 1 .3  cm.  The  size  of  the  plastic  zone  can  also  be  verified 
from  the  plot  of  isolines  of  the  <Tzz  component  of  stress  deviator  shown  in  figure  10. 

In  figure  1 1  the  calculated  trajectory  of  the  front  end  of  the  projectile  is  compared  with 
the  recorded  sequence  of  trigger  times  from  the  break  gages  in  the  study  in  reference  8. 
The  EOS  model  for  concrete  employed  is  shown  to  have  no  effect  on  the  calculated 
trajectory,  which  is  almost  identical  for  both  models.  The  trigger  times  from  the  gages  refer 
to  significant  deformations  exerted  on  the  copper  foils  that  are  sufficient  to  break  them 
causing  interruption  in  the  applied  voltage.  This  event  is  approximately  coincident  with  the 
time  of  arrival  of  the  elastic-plastic  boundary.  The  distance  between  the  regression  of 
trigger  times  and  the  trajectory  measures  the  size  of  plastic  zone  in  front  of  the  projectile, 
which  according  to  figure  1 1 ,  yields  approximately  Ap »  1 .3  cm.  As  seen  from  figure  9,  this 
measurement  of  the  plastic  zone,  correlates  well  only  with  results  obtained  using  the 
porous  EOS  model.  A  very  poor  agreement  in  the  size  of  the  plastic  zone  is  obtained  in 
case  of  the  analysis  using  the  sold  EOS  model.  Hence,  in  spite  of  relatively  "better" 
predictions  of  crater  size  with  the  solid  EOS  model  (both  calculations  significantly  disagree 
with  the  experimental  data  (fig.  5),  this  fact  provides  experimental  evidence  that  the  porous 
EOS  model  is  more  realistic  in  describing  the  effects  of  elastic-plastic  flow  in  concrete. 

The  initial  porosity  is  one  of  the  parameters  appearing  in  the  porous  EOS  model, 
and  the  sensitivity  of  penetration  calculation  to  the  values  of  are  presented  in  figures  12 
to  14.  The  measured  value  of  initial  porosity  of  Gregson's  (ref  14)  concrete  was  ^o=  0.18, 
and  this  value  was  adopted  as  a  baseline  in  this  work.  Figure  12  exhibits  hole  profiles  at 
500//S  calculated  with  values  of  initial  porosity  of  =  0.13,  <po  =0.18,  and  =  0.23,  and 

the  same  initial  density  po  =  2.2Ag/cm^.  The  values  of  other  parameters  for  the  concrete 
material  model  are  given  in  table  2.  Computational  experiments  indicated  that  varying  the 
value  of  initial  porosity  from  0.13  to  0.23  resulted  in  nearly  identical  penetration  depths  and 
hole  profiles.  Figure  13  presents  plots  of  isolines  of  porosity  calculated  for  varying  values 
of  (po,  which  suggest  a  cause  for  the  weak  correlation  between  the  calculated  hole  profiles 
and  the  initial  porosity.  As  projectiles  penetrate  the  target,  concrete  adjacent  to  the 
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boundary  of  the  forming  crater  is  crushed  and  all  porosity  is  squeezed  out.  Since  in  all 
cases  the  crushed  concrete  was  brought  to  the  same  state  with  ->  0,  the  further  motion 
of  this  material  is  nearly  identical,  which  results  in  similar  hole  profiles.  Physical  models  for 
concrete  can  be  approximated  by  a  composite  material  composed  of  a  crushable  cement 
matrix  and  an  aggregate.  Varying  the  value  of  initial  porosity  while  keeping  the  initial 
density  constant,  implies  materials  with  the  same  average  density  and  varying  size  of  the 
aggregate.  From  the  point  of  view  of  the  one-dimensional  hydrodynamic  theory  of 
penetration,  the  penetration  depth  is  a  function  of  the  target  material  density  and  since  the 
initial  density  was  the  same,  this  correlates  with  the  calculations  presented  in  figure  12. 

In  comparing  solid  and  porous  (^o  =  0.18)  EOS  models,  distinct  features  are 
observed  in  calculated  flow  variables  that  are  characteristic  to  each  of  these  models.  For 
the  numerical  validity  and  consistency  of  the  models,  it  would  be  reasonable  to  expect  that, 
with  decreases  in  the  initial  porosity,  the  characteristic  flow  variables  would  asymptotically 
approach  variables  calculated  with  the  solid  EOS.  As  representative  of  variables 
characteristic  for  the  flow,  radial  profiles  of  equivalent  plastic  strain  were  selected,  and  are 
shown  in  figure  14.  Comparison  of  these  plots  reveals  a  definite  tendency  for  gradual 
increases  of  the  size  of  the  plastic  flow  field  with  decreases  in  the  initial  porosity,  which 
transition  to  the  asymptotic  level  calculated  with  the  solid  EOS. 

Concrete  Targets:  Yield  Model 

Another  aspect  of  the  present  work  is  the  selection  of  a  material  model  for  concrete 
and  a  set  of  parameters  which  reproduce  the  experimentally  observed  penetration  depth 
and  the  profile  of  the  tunnel  portion  of  the  crater.  These  two  penetration  characteristics  are 
governed  primarily  by  a  constitutive  model  which  describes  the  elastic-plastic  behavior  of 
the  target  material  under  compression  exerted  by  the  projectile.  The  first  model  studied 
was  a  simple  model  of  an  ideal  elastic-perfectly-plastic  behavior  with  constant  yield 
strength.  This  model  has  only  one  adjustable  parameter,  (von  Mises  strength  Y)  which  can 
be  varied  to  match  the  experimental  data.  However,  it  is  not  physically  realistic  for  the 
description  of  the  complex  behavior  of  geologic  materials. 

Figure  15  shows  hole  profiles  calculated  with  a  constant  yield-strength  model  with 
varying  values  of  von  Mises  strength  Y.  The  unconfined  compressive  strength  of  the 
concrete  used  in  the  experiments  (ref  8)  was  0.374  Kbar,  and  in  the  numerical  experiments, 
the  value  of  von  Mises  yield-strength  was  varied  from  Y  =  0.1 87  Kbar  to  Y  =  1 .8  Kbar.  The 
results  represented  in  figure  15  indicate  that  the  calculated  hole  profiles  were  insensitive  to 
the  value  of  von  Mises  yield  in  the  constant  yield  model. 

The  validity  of  the  concrete  models  employed  in  the  numerical  analysis  presented  in 
this  work  is  established  through  comparison  with  the  following  experimental  data: 

•  The  size  of  the  elastic-plastic  zone  in  front  of  the  projectile 
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•  Hugoniot  data  from  the  shock  wave  studies  of  Gregson  (ref  1 4) 

•  The  resulting  profile  of  the  tunnel  portion  of  the  crater 

In  contrast  to  the  earlier  results  that  employed  a  constant  strength  yield  model,  figures  16 
and  1 7  present  the  results  of  calculations  based  on  a  model  with  a  pressure  dependent- 
yield  surface.  The  analysis  presented  in  the  earlier  section  showed  that  the  solid  EOS 
model  failed  to  reproduce  the  experimentally  observed  size  of  the  plastic  zone  in  front  of 
the  projectile.  This  result  is  in  agreement  with  the  shock  wave  studies  (ref  14),  which  show 
that  the  Hugoniot  data  for  concrete  cannot  be  reproduced  adequately  by  a  solid  EOS 
model.  Both  results  demonstrate  that  when  the  physical  mechanism  of  concrete  failure  due 
to  pore  closure  in  included  into  the  EOS  model,  only  then  the  structure  of  the  plastic  flow  is 
modeled  accurately.  Thus,  further  numerical  experimentation  with  a  solid  EOS  model  was 
abandoned.  All  further  calculations  employing  either  the  constant  yield-strength  model  or 
the  pressure  dependent  yield  model  were  performed  using  the  porous  EOS  model  only. 
Comparison  of  the  size  of  the  plastic  zone  in  front  of  the  projectile  calculated  in  all  these 
numerical  experiments  (table  2)  indicated  that  varying  the  value  of  the  von  Mises  yield- 
strength  in  the  constant-strength  model  did  not  affect  the  size  of  this  zone.  However,  using 
the  pressure-dependent  yield-strength  model,  not  only  was  the  experimental  depth  of 
penetration  reproduced,  but  also  the  measurements  of  the  profile  of  the  tunnel  portion  of 
the  crater  were  matched. 

The  configuration  of  the  yield  surface  Y  =  Y{p)  that  was  used  in  these  calculations 
was  adopted  from  Smith  (ref  16)  (figs.  18  and  19).  Figure  19  compares  the  calculated 
penetration  histories  of  V  =  const  and  Y  =  Y{p)  models.  According  to  the  one-dimensional 
hydrodynamic  theory  of  penetration,  the  penetration  into  the  target  is  achieved  through  the 
erosion  of  projectile  material.  With  the  higher  rates  of  erosion  of  projectile  material,  lower 
depths  of  penetration  are  expected.  The  velocities  of  the  front  end  of  the  projectile  are 
noticeably  different  for  the  two  models,  and  the  rate  of  projectile  erosion  for  Y  =  Y(p)  model 
is  greater  than  for  Y  =  const  model  (fig.  1 9). 


CONCLUSIONS 

The  structure  of  the  flow  field  around  the  penetrator  is  the  most  important  element 
that  governs  a  number  of  penetration  parameters  including  the  target  resistance  to 
penetration,  and  the  depth  and  the  diameter  of  the  produced  hole.  Since  concrete  is  a  very 
complex  material,  a  number  of  models  have  been  developed  to  describe  its  response  to 
high  rate  and  high  strain  deformations.  Since  the  choice  of  the  constitutive  model  affects 
the  elastic-plastic  zone  surrounding  the  penetrator,  several  candidate  models  were  studied 
using  three-dimensional  axisymmetric  analysis.  The  calculated  size  of  the  elastic-plastic 
zone  was  found  to  be  extremely  sensitive  to  the  choice  of  the  model.  Therefore,  the 
physical  validity  of  the  various  concrete  models  was  established  from  comparison  between 
the  analysis  and  the  existing  experimental  measurements  of  the  size  of  the  comminuted 
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concrete  zone.  From  all  the  models  tested,  only  the  porous  equation  of  state  model 
resulted  in  the  elastic-plastic  zone  size  predicted  by  the  experimental  dynamic 
measurements.  This  substantiated  the  theoretical  claims  that  the  porous  equation  of  state 
is  currently  the  most  realistic  model  for  representing  the  elastic-plastic  flow  in  the  concrete 
medium.  A  number  of  yield-strength  models  were  examined  for  modeling  the  constitutive 
behavior  of  concrete.  However,  the  pressure-dependent  yield  model  is  the  only  one  that 
resulted  in  hole  profiles  which  agreed  with  experiments.  Finally,  the  present  work 
established  a  reliable  material  model  for  concrete  targets  for  the  prediction  of  concrete 
penetration  for  a  wide  range  of  ballistic  impact  velocities. 


Spiral  reinforcing  bar 


Corrugated  steel  shell 


Figure  1 

Schematic  of  target  and  crater  resulting  from  projectile  impact 
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Figure  2 

Photograph  of  a  sectioned  target 
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(Spherical-nose  copper  projectile  impacting 
right  circular  cylinder  of  concrete  with  velocity 
Vo  =  0.1836  cm! /IS.) 


Figure  3 

Initial  configuration 
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Isolines  of  pressure  at  35  /js  for  solid  and  porous 


=  0.18)  EOS  models  for  concrete 
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Figure  5 

Hole  profiles  calculated  with  solid  and  porous 
{(fK)  =  0.18)  EOS  models  for  concrete 


50 
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Figure  6 

Profiles  of  axial  component  of  velocity  Vz  vs  axial  coordinate 
z  calculated  for  streamline  r  =  0  at  time  f  =  35  //s 


18 


Figure  7 

Profiles  of  radial  component  of  velocity  v,  vs  radial  coordinate 
r  calculated  at  z  =  5.0  cm  at  time  t=Z5  jus 
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(Calculated  at  time  t  =  35  jus) 


Figure  8 

Effect  of  EOS  model  for  concrete  on  configuration  of  flow  field 
in  front  of  projectile  -  isolines  of  equivalent  plastic  stream 
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Equivalent  plastic  strain 


(Calculated  at  time  t  =  35  ps) 


Figure  9 

Effect  of  EOS  model  for  concrete  on  configuration  of  flow  field  in  front 
of  projectile  -  profiles  of  equivalent  plastic  strain  vs 
axial  coordinate  z  for  streamline  r  =  0 
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(Calculated  at  time  t  =  35  jus) 


Figure  10 

Effect  of  EOS  model  for  concrete  on  configuration  of  flow  field  in  front 
of  projectile  -  isolines  of  component  of  stress  deviator 
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(The  distance  between  the  regression  of 
trigger  times  and  the  trajectory  of  the  front 
of  the  penetrator  measures  the  size  of  the 
plastic  flow  field  in  the  front  of  the  projectile.) 


Figure  1 1 

Trajectory  of  front  end  of  penetrator  compared 
with  trigger  times  from  break  gages 


23 


(Resulted  in  nearly  identical  penetration 
depths  and  hole  profiles.) 


Figure  12 

Analysis  with  porous  EOS  model  with  values  of  initial 
porosity  ^varying  between  0.13  and  0.23 
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(Concrete  adjacent  to  boundary  of  crater  was 
crushed  and  ail  porosity  was  squeezed  out.) 


Figure  1 3 

Isolines  of  porosity  at  time  t  =  35/xs  calculated  with 
different  values  of  initial  porosity  (Part  1  of  2) 
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(Concrete  adjacent  to  boundary  of  crater  was 
crushed  and  all  porosity  was  squeezed  out.) 


Figure  13 

Isolines  of  porosity  at  time  f  =  35  //s  calculated  with 
different  values  of  initial  porosity  (Part  2  of  2) 


Equtvaimt  plastic  strain 


r,  cm 


(Decreasing  the  value  of  initial  porosity  increases 
the  size  of  plastic  flow  field  in  front  of  projectile.) 


Figure  14 

Equivalent  plastic  strain  vs  radial  coordinate 
r  at  z  =  5.0  cm  at  time  /  =  35  //s 
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Figure  15 

Hole  profiles  calculated  with  constant  strength  model  for 
varying  values  of  von  Mises  strength  Y  (Part  2  of  2) 
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Pr^ssurt:  Log(p+3xi0‘*);  p,  Mbar 


Figure  18 

Pressure-dependent  yield  surface  Y  =  Y(p) 
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(Calculated  with  constant  strength  and 
pressure-dependent  strength  models.) 


Figure  19 

Trajectory  of  front  end  of  projectile 


Table  1 

Computations  performance  for  200  //s  of  penetration  time 


File  ID 

MeshMoUoa 

Mm/max  ibM  of  oonqMiUtkNMl  oellf ,  m 

No.  of 

Required 
CPU  time 

Zonel 

1  Zone  2 

Zone  3 

Zone  4 

compt. 

cells 

Arbax 

A**  in 
Ar.u 

A^in 

A^ 

A^bin 

AQik 

Ar«ax 

A^iin 

A^mx 

Ar*in 

Ar.ax 

hr 

pci0_00 

Euier  k  ALE 

0.0667 

2.1670 

0.0342 

p.1920 

0.0667 

0.6670 

0.0342 

1.7520 

0.2000 

0.2540 

0.0342 

1.7520 

9.2000 

9.2000 

0.0342 

1.7520 

29,895 

102.2 

pci0_01 

Euler  Si  ALE 

0.0650 

2.8250 

0.0722 

0.2000 

0.1333 

0.1333 

0.0722 

3.7190 

0.1400 

0.5080 

0.0722 

3.7190 

0.1333 

11.633 

0.0722 

3.7190 

10,971 

18.4 

pci0_02 

Euler  Si  ALE 

0.0650 

2.8250 

0.0722 

0.2000 

ITSsSQ 

0.2500 

JiM2 

3.7190 

0.5080 

0.0722 

3.7190 

0.2333 

11.633 

0.0722 

3.7190 

9,437 

9.4 

pci0_02E 

Euler 

r0.0650 

2.8250 

0.0722 

0.2000 

0.2500 

0.2500 

0.0722 

3.7190 

0.2330 

0.5080 

0.0722 

3.7190 

OM 

11.633 

0.0722 

3.7190 

9,437 

30.4 

(Benchmarked  on  a  SUN  4/490  workstation  for  varying 
degrees  of  discretization  and  modes  of  grid  motion.) 


Table  2 

Parameters  for  solid  and  porous  EOS  constitutive  models 


File  ID 

EOS  Form 

B 

a 

B 

B 

B 

TV) 

po 

Pr 

■a 

G 

Vmin 

pccu_02_8ol_fl 

SoUd 

Q 

.3122 

.4396 

.1693 

D 

D 

□ 

N/A 

2.240 

N/A 

0.374 

023 

0.9934 

Porous 

□ 

.3122 

.4396 

.1693 

□ 

D 

□ 

2.240 

2.732 

0.374 

0.9962 

pccu_02_por__f2 

Porous 

El 

.3122 

.4396 

.1693 

□ 

D 

□ 

0.23 

2.240 

2.909 

0.374 

0.9962 

pccu_02_por_13 

Porous 

Q 

.4396 

.1693 

□ 

D 

□ 

03 

2.240 

2.575 

0.374 

0.202 

0.9962 

pccu_03_Y_^f2 

Porous 

El 

.4396 

.1693 

a 

D 

El 

g[|3 

2.240 

2.732 

0.187 

0.202 

0.9969 

pccu_03_Y_f3 

Porous 

□ 

.4396 

.1693 

a 

D 

□ 

rim 

2.240 

2.732 

0.748 

0.202 

0.9969 

pccu_03_Y_f4 

Porous 

Q 

.4396 

.1693 

a 

D 

□ 

G&l 

2.240 

2.732 

1.496 

ggg 

0.9969 

pccu_03_Y_15 

Porous 

Qj 

.4396 

.1693 

D 

D 

□ 

123 

2.240 

2.732 

1.800 

0.9969 

Porous 

0 

.4396 

.1693 

D 

D 

□ 

nm 

2.240 

2.732 

Bia 

0.202 

0.9969 

(/Ai,  Si-1 ,  i  =  1 ,3,  and  G  are  Mbar,  Yo  is  in  Kbar,  po  and  pr  are  in  glcm^.) 
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